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ABSTRACT 

The  Effect  of  Pipe  Spacing  On  Marine  Pipeline  Scour 

(December  1988) 

Joseph  Henry  Westerhorstmann,  B.S.,  University  of  Florida 

Chair  of  Advisory  Committee:  Dr.  Jerry  L.  Machemehl 

Scour  is  a  major  threat  to  deep  water  pipelines.  A  study  was  performed  to 
determine  the  relationship  between  pipe  geometry  and  scour  depth,  and  examine 
scour  hole  dimensions  under  multiple  pipes.  Model  tests  were  conducted  in 
a  two-dimensional  wave-current  flume  facility  in  which  pipes  were  exposed  to 
unidirectional  flow  and  unidirectional  flow  with  oscillatory  motion.  The  model 
testing  investigated  pipes  in  contact  with  a  sand  bed  and  the  effect  of  pipe 
spacing  on  scour.  Scour  hole  profiles  were  measured  and  plotted  to  show  the 
scour  hole  shape  under  multiple  pipe  systems.  For  two  pipe  systems,  spacing 
between  pipes  of  1/2  pipe  diameter  resulted  in  less  scour  depth  than  full  diameter 
spacing.  Unidirectional  flow  with  oscillatory  motion  resulted  in  less  scour  than 
unidirectional  flow. 
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CHAPTER    I 
INTRODUCTION 


Purpose 

The  most  economical  method  of  transporting  oil  onshore  from  offshore 
production  facilities  is  through  submarine  pipelines.  Submarine  pipelines  are  also 
used  for  the  disposal  of  industrial  and  sanitary  wastes.  The  stability  of  these 
subsea  pipelines  is  threatened  by  loss  of  bed  support  resulting  from  scour.  This 
thesis  reviews  pipeline  scour  and  presents  results  of  model  testing  on  single  and 
multiple  pipes  in  contact  with  a  sand  bed. 

Few  theoretical  models  exist  which  accurately  describe  or  predict  scour  depth 
beneath  pipelines.  The  most  tvidely  accepted  models  predicting  scour  depth  are 
based  on  experimental  research. 

The  objectives  of  this  research  are  : 

(1)  To  determine  the  relationship  between  pipe  geometry  and  scour  depth, 
and 

(2)  To   examine   scour   hole   dimensions   under   multiple   pipes   subject    to 
undirectional  flow  and  unidirectional  flow  with  oscillatory  motion. 


The  citations  on  the  following  pages  follow   the  style  of  the   Journal  of 
Waterway,  Port,  Coastal  and  Ocean  Engineering  Division,  ASCE. 
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Pipe  spacing  and  scour  depth  were  studied  to  deternaine  the  difference  in  scour 
depth  between  a  single  pipe  and  multiple  pipes  subjected  to  similar  conditions. 
Single  and  multiple  pipe  systems  were  exposed  to  both  currents  and  currents  with 
wave  action  (oscillatory  motion). 

Problem 

Pipelines  are  being  placed  in  deeper  waters  over  greater  distances.  Pipelines 
on  the  ocean  floor  cause  disruptions  of  the  flow  pattern  and  changes  in  the  local 
fluid  velocity  which  create  bed  scour.  Damage  usually  occurs  during  severe  storms, 
when  large  ocean  waves  increase  bottom  currents.  Because  of  larger  horizontal 
drag  and  inertia!  forces  as  the  storm  intensifies,  scour  may  uncover  and  expose  a 
buried  pipeline.  The  exposed  pipe  may  be  damaged  since  pipelines  generally  are 
not  designed  for  conditions  where  spanning  and  vibrations  due  to  vortex  shedding 
occur  (Herbich  et  al.  1984). 

In  shaUow  depths,  pipelines  are  usually  buried,  placed  in  trenches  or  undergo 
self  burial  (Hulsbergen  1984).  A  severe  storm  generates  combined  current  and 
wave  action  on  the  seabed,  producing  higher  than  normal  fluid  velocities.  The 
combined  current  and  wave  action  effect  on  the  seabed  results  in  additional 
sediment  transport,  exposing  buried  pipelines.  Once  exposed,  the  pipe  lays  on 
the  surface  of  the  seabed,  and  is  further  subjected  to  scouring  currents. 


Local  and  area  scour  threaten  pipeline  stability  which  may  damage  the  pipe. 
Scour  holes  endanger  stability  when  large  free  spans  create  unacceptable  stresses 
due  to  static  deformations,  (Bruschi  et  al.  1986).  During  the  development  of  the 
scour  hole  under  the  pipeline,  the  stresses  inside  the  wall  of  the  pipe  increase  due 
to  the  sagging  of  the  pipe.  The  suspended  part  of  the  pipe  deflects  like  a  uniformly 
loaded  beam  w^ith  fixed  or  simply  supported  ends,  Fig.  1.  The  pipe  will  rupture 
if  the  static  bending  stresses  exceed  the  ultimate  strength  of  the  pipe. 

Flow  around  the  pipe  creates  a  series  of  complex  vortices  which  cause  the 
suspended  section  of  the  pipe  to  vibrate.  This  results  in  the  development  of  a 
deeper  scour  hole.  Pipe  vibrations  cause  pipe  failure  due  to  fatique.  Scour  is  a 
direct  consequence  of  vortex  shedding.  As  the  scour  hole  deepens,  the  frequency  of 
vibrations  may  increase  until  the  pipe  ruptures.  Vibrations  may  be  caused  by  the 
external  fluid  velocity  on  the  pipe  or  the  internal  fluid  transport,  or  a  combination 
of  the  two.  A  ruptured  pipeline  will  discharge  its  product  (such  as  oil)  into  the 
environment.  High  repair  costs,  income  loss  from  down  time  and  shorter  lifetime 
are  experienced. 
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CHAPTER    II 
LITERATURE  REVIEW 


Definition  of  Scour 

The  mechanism  of  scour  is  characterized  by  a  fine  balance  between  the 
local  flow  condition,  sediment  transportation  and  sediment  deposition.  Based 
on  research  conducted  by  Fredsoe  et  al.  (1988),  vortex  shedding  in  the  pipeline 
wake  is  the  the  major  cause  of  scour.  Turbulence  diffused  from  the  leeward  side  of 
the  pipe  wall  allows  sediment  to  remain  in  suspension,  allowing  for  easy  transport 
by  currents.  Scour  will  occur  if  the  rate  of  sediment  transport  upstream  of  the 
pipe  is  less  than  the  transport  downstream. 

Scour  depth  is  a  function  of  sediment  characteristics,  pipe  geometry,  flow 
and  wave  properties.  Among  all  the  variables  involved  in  the  phenomenon  of 
equilibrium  scour  depth,  scour  depth  (Sd)  is  the  dependent  variable  and  can  be 
expressed  as  a  function  of  the  following  independent  variables: 

Mean  fluid  velocity  {U), 

Critical  velocity  of  sediment  motion  (Vc), 

Water  depth  (y), 

Pipeline  diameter  (D), 

Wave  height  (//"), 


Wave  length  (X), 

Wave  period  (T), 

Dynamic  viscosity  of  water  (/i), 

Fluid  density  (p), 

Effective  sediment  size  {dg), 

Sediment  density  {ps), 

Relative  density  of  sediment  in  water  (A  =  (^^^^)),  and 

Gravitational  acceleration  (g). 

Note  that   kinematic  viscosity  {u   =   fi/p)  replaces  dynamic  viscosity  over 
density  of  water. 

Ibrahim  and  Nalluri  (1986)  developed  the  following  dimensionless  relation- 
ships between  variables: 


Sd/D  =  /[  {UD/u),  {Uy/u),  {U/y/^),  {H/gT'),  (UT/D),  (j) 

(HL'/y'),  (f7^/^),  (H/D),  {L/D),  {tj /V.),  {U / ./^Kd,)] 

The  most  significant  relationships  from  Eq.    1  which  influence  scour  are  the 
following  dimensionless  parameters: 


(UD/u)  and  {Uy/v)  are  pipe  and  flow  Reynolds  Numbers,  measuring 
turbulence, 

[H I gT^)  is  wave  steepness,  which  is  the  disturbing  capability  of  the  wave, 

{HL'^/y^)  is  the  wave  Ursell  Number  which  defines  the  relative  magni- 
tude of  wave  motion  and  water  depth. 

{U / \/gD),  {U j\^)  and  {U /\/gAds)  are  the  pipe,  flow,  and  sediment 
Froude  Numbers,  the  prime  force  behind  the  movement  of  the  seabed, 
and 

(UT/D)  is  the  Keulegan-Carpenter  number,  describing  the  dominant 
type  of  flow,  (i.e.,  current  or  wave). 


Incipient  Motion 

The  action  of  flow  on  the  bed  which  initiates  sediment  motion  is  measured 
by  the  critical  bed  shear  stress  {tc).  Critical  bed  shear  stress  is  expressed  as  the 
dimensionless  shear  stress  or  Shields  parameter  {6),  defined  as: 


(2) 


where: 

Tc  =  v;V  (3) 

The  critical  bed  shear  stress  Tc  is  related  to  the  bed  friction  by: 

Tc  =  cpU^o  (4) 

where: 

c  —  frictional  drag  coefficient  of  sand  bed,  and 
^100  —  fluid  velocity  at  100  cms  above  the  bed. 

Most  of  the  data  on  critical  shear  stress  comes  from  observations  made  in 
flume  experiments.  These  experiments  show  sediment  motion  as  highly  unsteady 
and  nonuniformly  distributed  over  the  bed  area.  In  near  critical  conditions,  grain 
motion  occurs  in  turbulent  bursts  over  small  areas,  increasing  as  near  shear  stress 
conditions  increase. 


When  shear  stress  is  near  the  critical  value,  gusts  of  sediment  motion  become 
random  in  space  and  time.  This  was  observed  by  Shields  (1936),  who  described  the 
initiation  of  motion  as  statistical  in  nature.  Einstein  (1942)  was  first  to  develop 
transport  relations  based  on  probabilistic  concepts.  Due  to  the  statistical  nature 
of  incipient  motion,  it  is  difficult  to  define  a  true  condition  for  initiation  of  motion. 

Critical  shear  stress  data  is  based  on  arbitrary  definitions  of  critical  condi- 
tions. \'anoni  (1966)  refers  to  the  following  three  intensities  of  motion  near  the 
critical  or  threshold  condition: 

(1)  Weak  movement —  few  or  several  of  the  smallest  sand  particles  are  in 
motion  in  isolated  spots  in  sufficiently  small  quantities  such  that  those 
moving  in  one  square  centimeter  of  the  bed  can  be  counted. 

(2)  •  Medium  movement —  grains  of  mean  diameter  are  in  motion  in  numbers 

too  large  to  be  counted.  The  movement  is  no  longer  local  in  character. 
It  is  not  sufficiently  strong  to  affect  the  bed  configuration  and  does  not 
resiilt  in  appreciable  sediment  discharge. 

(3)  General  movement —  sand  grains  up  to  and  including  the  largest  are  in 
motion  and  movement  is  occurring  in  all  parts  of  the  bed  at  all  times. 
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The  definition  of  critical  conditions  varies  widely  among  scientists  and  engi- 
neers. This  explains,  in  part,  the  variation  in  research  results  on  sediment  motion. 
However,  the  weak  movement  described  above  agrees  with  most  studies.  It  also 
agrees  with  Shields  critical  condition  (Shields  1936). 

As  the  velocity  of  fluid  flow  over  a  bed  of  sediments  is  increased,  the  fluid 
exerts  a  force  or  stress  on  the  particles.  Eventually  this  force  becomes  sufficient 
for  particle  movement  and  results  in  transporting  sediment  from  the  bed.  This 
is  known  as  the  threshold  of  sediment  movement,  or  the  critical  stage  of  erosion. 
The  fluid  velocity  which  results  in  threshold  movement  is  the  critical  velocity  of 
sediment  movement,  Vc- 

Threshold  curves  form  the  basis  for  predicting  initial  sediment  motion. 
Curves  based  on  a  combination  of  sediment  grain  diameter,  density,  velocity,  stress 
or  friction  are  entered  for  values  of  sediment  movement. 

Shields  diagram,  shown  in  Fig.  2,  is  a  relationship  between  the  boundary 
Reynolds  number,  (i2.)  and  Shields  parameter,  {6).  The  Shields  parameter  6  is 
defined  by  Eq.    2. 

The  boundary  Reynolds  number  is  defined  as: 


iJ-.  =  ^  ■  (5) 


where: 
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da  =    effective  sediment  diameter, 

u  —    kinematic  viscosity  of  water,  and 

U»  =    turbulent  shear  velocity. 

Turbulent  shear  velocity  [/»  is  defined  as: 

U.  =  ^J^p  (6) 

where: 

To  —  bed  shear  stress,  and 
p  —  density  of  water. 

Bed  shear  stress  r^  is  defined  as: 


Old] 


(7) 


where: 

dg  ■=  effective  sediment  size, 

Ci  =  form  coefficient  defining  the  effective  surface  area  of  the  particle, 

and 

Fd  =  drag  force  on  the  particle. 
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The  effective  surface  is  defined  as  the  projection  of  the  particle  on  a  plane 
perpendicular  to  the  direction  of  fluid  motion. 


Drag  force  on  a  particle  is  defined  as: 

Fd  -  ICopAUJ  (8) 

where: 

Co  —  drag  coefficient, 

-4  =  projected  area  of  particle  on  a  plane  perpendicular  to  the  fluid  flow 

direction, 

Ud  =  local  velocity  at  distance  d  above  bed,  and 

p  —  density  of  water. 


Shields  determined  this  relationship  by  measuring  bed  load  transport  for 
values  of  B.  Values  for  6  were  measured  at  greater  than  twice  the  critical  value.  He 
then  extrapolated  6  to  the  point  of  the  vanishing  bed.  This  procedure  was  used  to 
avoid  random  orientation  of  grains  and  variations  in  flow  conditions  which  result 
in  grain  movement  when  0  is  smaller,  Simons  and  §enturk  (1977). 

Shields  diagram  expresses  the  threshold  of  sedirhent  movement  as  a  function 
of  the  following: 


Sediment  density  (/Jj), 
Effective  grain  size  (f^,), 


Bed  shear  stress  {to). 
Critical  bed  shear  stress  (r^), 
Kinematic  fluid  viscosity  (i/),  and 
Fluid  density  (/>). 

The  determination  of  incipient  motion  is  important  in  the  study  of  pipeline 
scour.  Incipient  motion  of  sediment  leads  to  scour.  Shields  pubhshed  his  diagram 
for  critical  shear  stress  utilizing  boundary  layer  theory  and  fluid  mechanics. 
Considerable  dissatisfaction  with  this  shear  stress  criterion  is  presented  in  the 
Hterature,  Liu  (1958),  Vanoni  et  al.  (1966)  and  Yang  (1973). 

Yang  (1973)  summarizes  this  dissatisfaction  with  the  following: 

(1)  By  assuming  the  existence  of  a  universal  velocity  distribution  law, 
•  shear  velocity  or  shear  stress  is  a  measure  of  the  intensity  of  turbulent 

fluctuations.  Knowledge  of  turbulence  is  limited  mainly  to  model  tests. 

(2)  Shields  derived  his  criterion  on  incipient  motion  by  using  the  concept 
of  laminar  sublayer.  According  to  that  theory,  the  laminar  sublayer 
should  not  have  any  effect  on  the  velocity  distribution  when  the  shear 
Reynolds  number  is  greater  than  70.  Shields  diagram  indicates  that  the 
dimensionless  critical  shear  stress  still  varies  with  shear  velocity  Reynolds 
number  greater  than  70. 
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(3)  Shields  extends  his  curve  to  a  straight  line  when  the  shear  velocity 
Reynolds  number  is  less  than  3.  This  means  that  when  the  sediment 
particle  is  very  small,  the  critical  tractive  force  is  independent  of  the 
sediment  size.  It  has  been  shown  that  the  critical  tractive  force  is 
proportional  to  the  sediment  size. 

(4)  The  critical  shear  stress  cannot  be  determined  directly  from  Shields 
diagram.  It  must  be  determined  through  trial  and  error. 

(5)  Shields  simplifies  the  problem  by  neglecting  the  lift  force  and  considers 
tangential  force.  Lift  force  must  be  considered  for  high  shear  velocity 
Reynolds  number. 

(6)  Sediment  transport  cannot  be  uniquely  determined  by  shear  stress.  Using 
critical  shear  stress  as  the  criterion  for  incipient  motion  for  the  study  of 
sediment  transport  is  questionable. 

Based  on  the  preceding,  it  is  doubtful  that  the  Shields  diagram  is  the  best 
criterion  for  incipient  motion. 

Fig.  3  shows  curves  on  critical  velocity  plotted  against  mean  sediment  size 
for  quartz  sediment  in  water  {ps  =  2.65g/cm^)  as  presented  by  Vanoni  (1966). 
The  curves  shown  are  for  the  upper  and  lower  limit  of  the  mean  critical  velocity 
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taken  from  several  sources.  The  curves  are  for  depths  of  at  least  one  meter.  The 
relationship  was  developed  by  fitting  the  curves  with  observed  data.  The  bottom 
velocities  are  obtained  by  extrapolating  velocity  measurements  to  the  plane  of  the 
bed. 


Scour  Depth  and  Profile 

Kjeldsen  (1973)  proposed  the  following  formula  for  scour  depth: 

5(i -0.972(f7Vp)°'  D'-'  (9) 

where: 

Sd  =  maximum  equilibrium  scour  depth, 
('  =  mean  local  fluid  velocity, 
g  =  gravity,  and 
.  D  =  pipe  diameter. 

Bijker  and  Leeuwestein  (1983)  proposed  the  following  modification  to  Kjeld- 
sen's  equation,  w^hich  incorporates  mean  grain  size: 

Sd  =  0.929(C^^/25)°"  D'-''  (d5o)"°°'  (10) 

where:     dso  =  mean  particle  diameter. 
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Mean  particle  diameter  is  the  grain  diameter  of  the  sediment  particle  such 
that  half  the  sediment  particles  (by  weight)  have  larger  diameters  and  half  have 
smaller  diameters. 


Mao  (1986)  developed  the  scour  hole  profile  from  potential  flow  theory  as: 


fe+f 


+  ^{l+cos(2jri))     for  III  < 
for  \x\  > 


(11) 


where: 

h{x)  =  Vertical  distance  from  pipe  center  to  bottom  of  scour  hole, 

w  =  Width  of  trench,  and 

€  =  Distance  between  pipe  and  bed. 
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Experimental  Results 

The  most  recent  significant  model  test  on  pipe  scour  is  based  on  the  inter- 
action between  a  fixed  pipeline  and  an  erodible  bed.  Mao  (1986)  studied  the 
influence  of  a  fixed  pipe  on  an  erodible  sand  bed  under  various  wave  and  current 
conditions.  The  experimental  set-up  was  typical  for  model  tests  on  pipes.  Tests 
were  run  on  both  a  5  cm  pipe  and  a  10  cm  pipe  placed  on  a  quartz  sand  bed. 

The  model  test  resulted  in  the  following  conclusions: 

(1)  The  scour  process  consists  of  a  jet  period  and  a  wake  period, 

(2)  The  flow  stagnation  point  is  on  the  upper  part  of  the  pipe, 

(3)  The  scour  rate  is  very  high  at  an  early  stage,  later  on  it  decreases, 

(4)  In  equilibrium,  bed  velocities  are  similar  to  mean  current  velocity, 

(5)  The  bigger  the  initial  gap,  the  shallower  the  scour  depth, 

(6)  The  scour  depth  is  a  weak  function  of  Shields  parameter,  0, 
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(7)      The  final  scour  depth  Sd  is  usually  less  than  the  diameter  of  the  pipe, 
and 


(8)      The  final  scour  depth  is  influenced  by  the  wake  behind  the  pipe. 

Fredsoe  et  al.  (1988)  studied  the  interaction  between  a  vibrating  pipe  and 
erodible  bed.  The  experimental  set-up  was  similar  to  the  preceding  model  test, 
with  the  following  difference:  a  spring  supported  rigid  cyhnder  was  used  as  the 
model  pipe,  which  was  then  vibrated. 

The  influence  of  pipe  vibration  on  scour  and  the  influence  of  scour  on  pipe 
vibration  results  in  the  following  conclusions: 

(1)  Pipe  vibration  induces  additional  erosion,  which  gives  rise  to  relatively 
larger  scour  depths  and  scour  widths.  This  effect  was  found  to  be  most 
pronounced  for  smaller  gaps  between  the  pipe  and  the  bed. 

(2)  Vibrations  are  dominated  by  vortex  shedding  due  to  this  extra  erosion. 
Additional  erosion  occurs  even  if  the  pipe  is  placed  very  close  to  the 
original  undisturbed  bed. 

Ibrahim  and  Nalluri  (1986)  studied  pipeline  scour  due  to  currents  and  currents 
with  waves.    Tests  were  performed  on  three  diff"erent  grades  of  bed  sand  of  sizes 
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dj  —  0.425  mm,  0.80  mm,  and  1.50  mm.  Pipe  sizes  were  2.50  cm,  5.00  cm,  and 
7.50  cm.  The  pipes  were  fully  exposed  to  flow,  i.e.  resting  on  a  sand  bed.  They 
also  studied  pipelines  that  were  one-half  and  three-quarters  buried.  The  action 
between  pipeline,  seabed  and  flow  due  to  the  combined  wave  and  current  was  more 
intense  compared  to  that  due  to  unidirectional  current. 

Ibrahim  and  Nalluri  (1986)  developed  dimensionless  plots  between  the  relative 
scour  depth,  Sd^D  and  the  following  dimensionless  parameters: 

(1)  [V !\fgD)  —  pipe  Froude  Number,  {Nfp),  Fig.  4, 

(2)  (r/v'^Affa)  —  sediment  Froude  Number,  (iVs)),  Fig.  5, 

(3)  {HL^ /y^)  —  wave  Ursell  Number,  Fig.  6,    and 


(4)     {H/gT'^)  —  wave  steepness,  Fig.  7. 


Figs.  4-7  attempt  to  identify  the  influence  of  several  groupings  in  Eq.  1  on 
relative  scour  depth,  Sd/D.  These  plots  become  useful  in  comparing  the  data 
obtained  in  this  research  to  that  of  Ibrahim  and  Nalluri  (1986). 
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Fig.  4.  Influence  of  pipe  Froude  Number   {tj fy/gB)  on  scour. 
(Ibrahim  and  Nalluri  1986) 
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Fig.  6.  Influence  of  wave  Ursell  Number  {HL^/y^)  on  scour. 
(Ibrahim  and  Nalluri  1986) 
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Despite  the  scatter  of  data  points  in  Fig.  4  and  5,  a  trend  reveals  a 
relationship  between  SdjD  and  flow  velocity.  The  flow  velocity  is  part  of  the 
Froude  Number,  which  represents  an  inertia-gravity  force  ratio.  The  relative 
scour  depth  Sd/D  increases  with  increasing  Froude  Number.  This  is  expected  — 
larger  fluid  velocity  creates  deeper  scour. 

Fig.   6  is  the  influence  of  the  wave  Ursell  Number  [HL"^ /y^)  on  scour.   The 
Ursell  Number  represents  relative  magnitude  between  H I L  and  y/i,  where: 
H  —  wave  height, 
L  =  wave  length,  and 
y  =  water  depth. 

The  region  where  f>  <  1  corresponds  to  deeper  water  and  a  region  where  Stokes 
wave  theory  is  apphcable.  This  group  tends  to  reveal  a  non-interdependency  of 
data.  According  to  Ibrahim  and  Nalluri  (1986),  the  wave  Ursell  Number  defines 
the  magnitude  of  energy  generated  by  the  wave  and  transmitted  to  the  bed.  The 
data  in  Fig.  6  fall  in  a  region  where  U^  >  ^t  signifying  shallow  water. 

Fig.  7,  the  wave  steepness,  describes  the  disturbing  capability  of  the  wave 
on  the  sediment.  For  shallow  depths,  the  data  indicates  an  increase  in  Sd/D  for 
steeper  waves. 
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CHAPTER    III 
PIPELINE  SCOUR 


Types  of  Pipeline  Scour 

Two  categories  of  scour  are  global  and  local.  Global  scour  is  an  area  erosion 
of  the  bed  in  the  vicinity  of  the  disturbance.  Generally,  it  is  an  eUiptical  shape 
with  the  major  axis  in  the  direction  of  the  dominant  flow.  Local  scour  develops 
at  the  seabed  immediately  adjacent  to  the  pipeline.  Local  pipeline  scour  results 
in  a  scour  hole  below  the  pipeline  (Dahlberg  1983). 

When  examining  local  scour,  the  two-dimensional  case  is  considered.  The 
primary  cross-sectional  area  is  a  transverse  cut  through  the  pipe  and  seabed. 
Span  analysis  examines  the  cross-sectional  area  along  the  pipe  axial  length.  Flow 
parameters  which  influence  scour  are  steady  flow  (current),  oscillatory  motion 
(wave),  or  a  combination  of  current  with  superimposed  wave  velocity.  The  net 
effect  is  a  pulsating  mass  transfer  in  one  direction.  The  combined  oscillatory 
motion  and  current  velocities  yield  the  mean  flow  in  one  direction.  The  scour 
due  to  a  steady  current  is  greater  than  the  combined  wave-current  affect. 


26 

Leeuwestein  (1984)  describes  three  types  of  local  scour  based  on  relative 
position  to  the  pipe  as  tunnel,  luff  and  lee  erosion,  Fig.  8.  Lee  erosion  occurs  while 
the  sediment  remains  in  suspension  in  the  wake.  During  lee  erosion  the  sediments 
are  transported  up  by  jet  flow  into  the  wake  and  then  carried  away  by  currents. 
Total  scour  depth  is  larger  for  current -induced  scour  than  for  wave-induced  scour. 
Bijker  and  Leeuwestein  (1983)  found  that  wave  induced  scour  is  about  30%  the 
depth  of  current-induced  scour.  The  reason  for  smaller  wave-induced  scour  is 
that  the  orbital  nature  of  the  flow  results  in  some  sediment  remaining  in  the  scour 
hole. 

Lee  erosion  extends  further  downstream,  resulting  in  a  wide  trench  around 
the  pipe,  Fig.  8.  Luff  erosion  generally  appears  just  upstream  of  pipelines  in 
the  absence  of  tunnel  erosion,  usually  near  partially  buried  pipelines.  Luff  erosion 
is  not  of  major  concern  because  it  generally  does  not  result  in  free  spans.  The 
primary  influence  on  scour  beneath  a  pipe  is  lee  erosion.  Tunnel  erosion  occurs 
in  the  early  stages  of  the  scour  process,  then  lee  erosion  takes  over  during  the 
growth  phase  of  the  scour  hole.  Tunnel  erosion  is  the  scour  which  occurs  directly 
beneath  the  pipe.  Tunnel  erosion  leads  to  lee  erosion  which  results  in  free  spans. 
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Contributing  Factors 

Maximum  equilibrium  scour  depth  is  a  function  of  several  factors.  The  major 
cause  of  erosion  is  the  vortex  shedding  in  the  wake.  Vortex  shedding  is  turbulence 
diffused  from  the  leeward  side  of  the  pipe  wall,  allowing  sediment  to  remain  in 
suspension.  Sediment  is  transported  up  by  jet  flow  into  the  wake  of  the  vortex 
and  then  carried  off  by  currents.  Vortex  shedding  drives  lee  erosion. 

The  Keulegan-Carpenter  number  {KC)  is  the  parameter  which  describes  the 
type  of  flow,  being  either  steady  flow  (current),  oscillatory  motion  (wave),  or  a 
combination  of  the  two.  The  KC  number  is  a  non-dimensional  measure  of  the 
flow,  described  earHer  as: 

KC  =  UT/D  (12) 

where: 

U  =  local  fluid  velocity  near  pipe, 
T  =  wave  period,  and 
D  =  pipe  diameter. 

A  low  KC  number  (i.e.  <  50)  means  that  oscillatory  flow  dominates.  As  KC 
increases,  current  flow  dominates. 
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According  to  Fredsoe  et  al.(1988),  KC  represents  the  effect  of  the  lee-wake  on 
the  scour  and  constitutes  the  most  important  parameter  which  governs  equilibrium 


Mechanics  of  Scour 

Sediment  transport  conditions  do  not  necessarily  exist  upstream  of  the  pipe. 
The  presence  of  the  pipe  is  responsible  for  generating  enhanced  flow  conditions 
that  trigger  the  scouring  process.  A  pipeline  near  an  erodible  bed  is  subjected  to 
constantly  changing  scour  conditions.  Pipe  movement  above  the  scour  hole,  from 
sag  and  vortex  shedding  induced  vibrations,  also  contributes  to  the  variability  of 
flow  regimes. 

Sediment  mechanics,  incipient  motion  and  pipeline/soil  interaction  are  pri- 
mary factors  that  influence  scour.  Sediment  motion  is  determined  by  the  bed 
shear  stress,  which  is  defined  as  the  tangential  force  exerted  along  the  seabed  by 
the  passing  flow  of  water.  Scour  is  effected  by  the  leeward  wake  created  by  vortex 
shedding  of  the  pipe.  Hydrodynamic  forces  on  the  pipeline  may  influence  the 
scour  beneath  a  pipeline,  but  they  are  not  the  principle  driving  force.  Lift  forces 
or  pipe  vibrations  on  the  pipeline  raise  the  pipe  sufficiently  for  jet  flow  through  a 
gap  between  the  pipe  and  bed.  This  is  known  as  the  nucleation  phase  or  the  jet 
period.  During  this  phase  a  small  hole  is  created  under  the  pipe. 
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Incipient  motion  occurs  when  the  jet  velocity  between  the  pipe  and  bed 
exceeds  the  critical  velocity  needed  to  begin  sediment  movement,  and  tunnel 
erosion  is  initiated,  Fig.  9.  Vortex  shedding  creates  a  turbulent  wake  leeward 
of  the  pipe,  which  keeps  the  sediment  suspended  for  sufficient  time  for  current  to 
carry  the  sediment  downstream.  The  growth  phase  or  wake  period  of  the  scour 
hole  occurs  while  the  vortex  shedding  creates  a  turbulent  wake  leeward  of  the 
pipe. 

The  sediment  in  suspension  is  balanced  against  gravitational  forces  that  tend 
to  pull  it  back  into  the  bed.  Lee  erosion  takes  place  while  the  sediment  is  carried 
off  in  the  wake  of  vortex  shedding.  Vortex  shedding  may  be  further  enhanced  by 
pipe  vibrations. 

Equilibrium  scour  depth  is  reached  when  the  velocity  near  the  bed  surface 
is  about  equal  to  the  current  velocity.  Sediment  removal  from  the  scour  hole 
continues  as  long  as  the  critical  velocity  of  sediment  movement  exists.  Critical 
velocity  dominates  as  long  as  the  the  current  velocity  is  greater  than  the  threshold 
value  of  sediment  movement.  The  growing  scour  hole  alters  the  flow  until  a  new 
dynamic  equilibrium  is  reached  between  the  erosion  and  deposition.  The  increase 
in  the  scour  hole  size  causes  a  decrease  in  the  jet-like  effect  until  it  returns  to 
normal  mean  fluid  velocity.  This  happens  when  the  depth  of  the  scour  hole  is 
approximately  one  pipe  diameter  deep.  At  this  point  the  scour  process  ceases. 
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The  presence  of  the  pipeHne  causes  vortex  shedding.  Vortex  shedding 
complicates  evaluation  of  boundary  conditions  for  scour  initiation.  Boundary 
conditions  are  determined  by  the  type  of  soil.  Boundary  layer  parameters  are 
the  bed  shear  stress  (tq)  and  bed  friction  coefficient  (c),  both  parameters  of  the 
sediment.  Shear  stress  (r)  due  to  current  flow  or  oscillatory  motion  is  determined 
as  a  function  of  bed  roughness  and  effective  grain  size.  Critical  shear  stress  (tc) 
which  contributes  to  initial  movement,  is  a  function  of  the  immersed  weight  of  a 
grain  of  sediment,  mean  grain  size  and  critical  velocity. 

Time  scale  effect  on  the  formulation  of  scour  rate  is  related  to  the  dominant 
type  of  flow  through  the  Shields  parameter  {9).  The  scour  rate  is  high  at  first,  then 
it  decreases  with  time.  The  larger  the  Shields  parameter  (9),  the  higher  the  scour 
rate  in  the  early  stages.  The  presence  of  superimposed  waves  greatly  increases  the 
rate  at  which  sediment  is  transported,  which  accelerates  the  formation  rate  of  the 
scour  hole.  However,  the  equilibrium  scour  depth  is  less  in  the  combination  flow 
condition  than  in  steady  flow. 

The  shape  of  the  scour  hole  will  vary  in  symmetry,  and  is  also  a  function  of 
the  dominating  type  of  flow.  This  is  best  described  by  the  Keulegan-Carpenter 
( A'C)  number.  The  width  of  the  hole  is  a  function  of  the  sediment  characteristics, 
and  the  symmetry  is  a  function  of  the  type  of  flow.  In  a  unidirectional  flow,  the 
leeward  side  of  the  pipe  will  have  less  of  a  slope  than  the  upstream  side.  In  a 
two-directional  oscillatory  flow,  the  profile  is  symmetrical.  Both  the  symmetry 
and  dominating  flow  are  closely  related  to  the  Keulegan-Carpenter  {KC)  number. 
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Scale  Effects 

One  of  the  difficulties  in  conducting  model  tests  is  insuring  that  scale  factors 
are  accurate.  Scale  effects  are  based  on  scaling  laws  for  local  scour  and  are  related 
to  the  similarity  of  sediment  size,  grain  shape,  and  grain  Froude  number.  The 
limiting  grain  size  for  permissible  reduction  of  sediment  in  a  model  is  about  0.05 
mm.  Pipe  sizes  are  a  function  of  the  test  flume  dimensions.  Typical  offshore  oil 
pipelines  are  one  meter  in  diameter.  In  this  case  a  model  pipe  of  ten  centimeters 
represents  a  1:10  scale. 

Test  conditions  which  govern  the  appUcability  of  experiments  are  the 
Reynolds  [Re]  number  and  Keulegan-Carpenter  [KC]  number.  The  R^  num- 
ber is  a  measure  of  the  turbulence  and  the  KC  number  describes  the  type  of  flow. 
A  high  /?e  puts  the  flow  in  the  turbulent  range,  representing  typical  field  condi- 
tions. A  KC  number  greater  than  50  describes  a  steady  flow,  whereas  a  lov/er 
value  takes  into  account  the  oscillatory  flow  contribution  from  the  wave  velocity. 
For  wave-current  experiments  the  following  R^  and  KC  ranges  are  recommended: 


(1)     Reynolds  number: 


fle  =  ^>10^ 


(2)     Keulegan-Carpenter  number: 


UT 
KC=  —  <  50 
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CHAPTER    IV 
COMPUTER  PROGRAM 

Program  Development 

The  program  SCOUR. FOR  is  a  FORTRAN  coded  program  which  calculates 
the  maximum  scour  depth  and  scour  profile  based  on  equations  defined  earlier: 

(1)      Equilibrium  scour  depth,  (see  Fig.  10): 


where:    A',  is  a  sediment  factor  taking  into  account  effective  sediment  size 
{dg)  and  densicy  (^5). 


(2)     Scour  profile,  Eq.    7,    (see  Fig.  11): 


fe+f- 


,,    .       /  e  +  f  +  f  (1  +  co*(27ri))       for  |:.i  <  f 
''(^)  =  1  -  .   o  for  |x|  >  f 
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Fig.  10.  Definition  sketch  of  scour  under  a  pipe. 


W- 


Fig.  11.  Scour  hole  profile,  from  Mao  (1986). 
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SCOUR. FOR  inputs  Stokes  III  wave  theory  as  a  subroutine,  WAVE123. 
Stokes  V  wave  theory  is  the  recommended  guideline  adopted  by  the  American 
Petroleum  Institute  (API)  for  use  in  offshore  design  problems.  The  theory  has  wide 
appHcations  in  the  coastal  and  offshore  industry.  For  purposes  of  this  program, 
Stokes  III  will  suffice. 


The  WAVE123  subroutine,  modified  for  the  SCOUR. FOR  program,  computes 
the  maximum  wave  velocity  at  a  depth  of  interest  near  the  pipe.  The  wave  velocity 
is  combined  with  current,  yielding  the  mean  local  fluid  velocity  {U)  acting  on  the 
pipe.  The  input  wave  height  and  period  are  for  extreme  events,  resulting  in  the 
highest  bottom  velocities.  Extreme  events  are  storm  or  hurricane  conditions. 

The  program  is  governed  by  the  following  assumptions: 

(1)  The  maximum  wave  velocity  at  the  depth  of  interest  is  algebraically 
summed  with  the  current  velocity,  yielding  mean  fluid  velocity  acting  on 
the  pipe. 

(2)  Although  wave  velocity  is  oscillatory,  the  net  effect  is  a  pulsating  mass 
transfer  in  one  direction. 

(3)  A  normal  angle  of  attack  on  the  pipe  has  the  greatest  effect  on  scour 
depth,  resulting  in  maximum  scour  depth. 
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(4)        The   computer   program  is  based   on   a  unidirectional  flow   from  the 

combined   effect    of  current    and   wave   velocities.      The   velocities   are 

algebraically  added.     The  shape  factor  is  determined  by  the  width  of 

the  scour  hole  in  terms  of  pipe  diameters. 


(5)      The  program  is  designed  for  depths  in  the  100  meter  range  under  extreme 
storm  conditions. 


Program  Flexibility 

The  SCOUR. FOR  program  flow  chart,  Fig.  12,  is  a  visual  representation 
of  the  solution  sequence.  The  flow  chart  depicts  the  flow  of  control,  identifying 
interrelationships  and  decision  requirements  of  the  program  sequence.  The  flow 
chart  is  adjustable  to  other  computer  languages  and  systems  capable  of  executing 
basic  nurnerical  computations. 

The  program  accommodates  any  type  of  pipe  orientation  to  the  bed.  The 
gap  distance  value  e  must  be  input  if  the  pipe  does  not  rest  on  the  seabed.  When 
the  pipe  is  in  contact  with  the  seabed,  e  =  0.  The  scour  profile  coordinate  system 
is  measured  from  the  center  of  the  pipe  if  e  >  0,  otherwise  the  coordinate  system 
is  adjusted  to  the  seabed  for  €  =  0. 
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SCOUR. ROR 
CRI.L.  MRVEX23 
F^OR  Uui 

RERO 

1 

^ 

WRVE123 
CSTOKES  ZZZ3 
Uui  Q  Z 

RERD 

Fig.  12.  SCOUR.FOR  flow  chart. 
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The  effect  of  sediment  characteristics  is  controlled  by  the  sand  factor  (5/) 
and  the  sediment  factor  (Kg).  The  sediment  factor  will  influence  the  equiHbrium 
scour  depth  and  is  a  function  of  the  type  of  sediment.  It  is  adjustable  in  the 
program  listing. 

The  inputed  sand  factor  {Sf)  relates  the  pipe  diameter  to  the  scour  hole  width 
(u-),  and  is  a  function  of  Shields  parameter  {9).  The  symmetry  of  the  scour  hole 
is  controlled  in  the  program  list  and  is  a  function  of  the  dominant  type  of  flow, 
being  either  oscillatory  or  unidirectional.  The  Keulegan-Carpenter  (KC)  number 
governs  the  type  of  flow.  The  SCOUR.FOR  FORTRAN  program  CODE  Hst  is 
included  in    Appendix  I. 

Execution 

SCOUR.FOR  reads  pipe  diameter  (D),  sand  factor  (Ks)  and  current  velocity 
(f'c),  ^n  metric  units,  from  the  input  file  MAIN. DAT.  The  program  subroutine 
WAVE123  input  is  wave  height  (H),  period  (T),  water  depth  [d)  and  depth  of 
interest  near  the  pipe  {Z). 

SCOUR.FOR  program  output  creates  an  output  file  PROFILE.DAT,  which 
is  easily  plotted  using  any  standard  graphics  package. 
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The  program  was  designed  on  a  Digital  8000  VAX/VMS  mainframe  com- 
puting system  with  FORTRAN  compiler  located  at  Texas  A&M  University.  The 
SCOUR. FOR  and  WAVE123  FORTRAN  codes  are  loaded  into  the  mainframe 
and  executed  with  the  following  commands: 

(1)  FOR  SCOUR,WAVE123 

(2)  LINK  SCOUR,WAVE123 

(3)  RUN  SCOUR 

The  program  is  best  demonstrated  by  varying  the  MAIN. DAT  input  data 
(current  velocity  (Tc),  pipe  diameter  {D)  and  sand  factor  (5/)),  while  keeping  the 
\VAVE123  data  constant  (except  for  depth  of  interest  (y)  when  changing  (-D)),  and 
comparing  the  plotted  results.  The  sample  output,  Fig.  13,  was  created  from  the 
SCOUR. FOR  output  file,  PROFILE.DAT,  using  the  VAX/VMS  graphics  system 
PICSURE  package. 
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CHAPTER  V 
EXPERIMENTAL  APPARATUS  AND  EQUIPMENT 


Model  testing  was  conducted  in  a  two-dimensional  wave  and  flume  facility  as 
shown  in  Fig.  14.  The  wave-flume  is  45.72  m  long,  0.46  m  wide,  and  1.22  m  deep. 
One  side  of  the  flume  consists  of  a  plexiglas  side  wall  for  observing  tests.  Fig.  15 
is  the  schematic  layout  of  the  test  facility. 

A  6.0  m  long  ramp  section  consisted  of  a  3.0  m  long  gradual  rise  (1:15) 
followed  by  a  level  3.0  m  long  section,  Fig.  16.  The  ramp  preceded  the  sand  test 
section  and  provided  for  a  smooth  transition  from  the  flume  bottom  to  sand  bed. 
It  also  contributed  to  the  development  of  the  boundary  layer. 

A  test  section  located  approximately  20.0  m  from  the  wave  generator  con- 
sisted of  a  6.0  m  long  by  0.20  m  deep  sand  bed.  This  test  section  was  12  m  long, 
including  the  ramp  and  sand  bed  section,  shown  in  Fig.  16.  Pipes  were  placed  on 
the  bed  at  about  3.0  m  downstream  of  the  ramp.  This  allowed  for  development 
of  a  boundary  layer  and  equilibrium  condition  in  the  sand  bed  section  before  the 
current  reached  the  pipes. 
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Fig  14.  Test  facility. 
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Currents  were  produced  by  a  pumping  system  capable  of  developing  a  flow 
velocity  of  1.20  m/sec/  The  flow  rate  is  controlled  by  two  flow  control  valves. 
Water  depth  is  controlled  by  an  adjustable  weir  located  at  the  end  of  the  tank. 
Currents  were  controlled  by  adjusting  the  flow  control  valves  and  weir. 

The  wave  generator  produces  a  wave  by  the  plunging  motion  of  a  wedge, 
where  the  height  of  the  wave  is  determined  by  the  depth  of  the  plunger  and  the 
depth  of  water.  Wave  period  and  wave  length  were  adjusted  by  a  rheostat  which 
controls  the  plunger  speed.  The  wave  generator  was  preset  for  test  conditions.  A 
wire  mesh  wave  filter  which  dissipated  the  smaU  surface  ripples  from  the  generated 
waves  was  located  2.0  m  from  the  plunger. 

A  wave  absorber  was  placed  at  the  end  of  the  tank  to  reduce  wave  reflection. 
The  wave  absorber  design  permits  passage  of  large  volumes  of  water.  This  allows 
for  significant  currents  and  waves  with  undisturbed  current  flow  and  reduced  wave 
reflection. 

Three  sizes  of  pipes  are  utilized:  3.0,  6.0,  and  9.0  cm.  The  pipes  were  secured 
in  the  tank  through  compression  springs  and  dowels  placed  inside  each  pipe.  Each 
pipe  was  placed  in  the  tank  under  flow  conditions,  and  set  in  place  with  the  flow 
direction  normal  to  the  longitudinal  axis  of  the  pipe.  Proper  pipe  spacing  was 
achieved  through  the  use  of  spacers.  The  pipe  spacing  is  the  clear  width  between 
outer  pipe  diameters. 
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Current  velocity  was  measured  with  a  propeller  current  meter.  Velocity 
raeasurements  were  taken  at  5.0  cm  from  the  bottom  up  to  the  water  surface. 
Velocity  profiles  were  developed  by  measuring  the  current  at  5.  10,  20,  30,  40,  50, 
and  60  cm  from  the  bottom. 


A  Seasim  resistance  type  autocompensating  wave  gauge  and  Hewlett  Packard 
Model  17403 A  strip  chart  recorder  were  used  to  measure  and  record  the  wave 
characteristics.  The  wave  gauge  consists  of  a  40.0  cm  long  stainless  steel  probe 
with  anodized  aluminum  calibration  rod  and  clamp  attachment.  The  maximum 
recommended  cable  and  circuit  resistance  is  1.0  ohm  and  frequency  range  in 
autocompensation  mode  is  0.2  to  above  10.0  Hz. 

Scour  profiles  were  measured  with  a  point  gauge  profiler.  The  profiler  was 
manually  operated  with  a  precision  of  5.0  mm. 
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CHAPTER  VI 
EXPERIMENTAL  PROCEDURES 


Seventy-two  test  runs  were  conducted  for  this  study.  Pipe  sizes  of  3.0,  6.0, 
and  9.0  cm  were  placed  on  a  sand  bed  test  section.  Fig.  17  is  the  result  of  a 
sieve  analysis  of  the  sand;  it  shows  that  dso  —  0.08  mm.  Table  I  shows  the  test 
pipe  configurations,  flow  conditions  and  test  numbers.  Each  pipe  configuration 
was  subjected  to  four  flow  conditions: 


(1)     Current  only: 


(2)      Current  with  wave: 


U  =  0.37  m/sec 


U  —  0.37  ml  sec 


wave  characteristics: 

H  =  6.0  cm,    I  r=  1.0  m,    T  -  1.2  sec 
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Table  I.  Test  Numbers  and  Test  Conditions. 

TEST 

r   NUMBERS 

TEST 
CONnCURATlON 

U  m  0.37  m/«      U  m  0.2B  m/« 

PIPE  DIAMETER  -  0  -  (cm) 

6 

9 

3 

6 

9 

—  O 

1 

15 

29 

43 

53 

63 

—  o 

2 

16 

30 

44 

54 

64 

—  OIO 

3 

17 

31 

45 

55 

65 

—  OfO 

4 

18 

32 

46 

56 

66 

—  o»o 

5 

19 

33 

47 

57 

67 

—  O'O 

6 

20 

34 

48 

58 

68 

—  o?oio 

7 

21 

35 

49 

59 

69 

OlOfO 

8 

22 

36 

50 

60 

70 

—  0»  O"  o 

9 

23 

37 

51 

61 

71 

^=o  »o»  o 

10 

24 

38 

52 

62 

72 

—  UfCJ 

o 

11 

25 

39 

NOTES  ! 

—  0 

currant  only 

o 

12 

26 

40 

—  0»0 
n 

13 

27 

41 

currant  It  wav« 

U 

28 

42 

h«]ght  m  §  em 
langth  «  1.0  m 
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(3)      Current  only: 


(4)      Current  with  wave: 


U  =  0.26  m/sec 


U  =  0.26  m/sec 


wave  characteristics: 

/T  =  6.0  cm,    L  =  1.0m,    T  ^  1.2  sec 

where:  U  —  mean  local  velocity  5.0  cm  above  bed. 

The  pipe  spacing  is  the  clear  width  between  outer  pipe  diameters.  The  flow 
control  valve,  adjustable  weir  and  wave  maker  were  calibrated  and  preset  for  the 
test  conditions. 

The  critical  velocity  of  incipient  motion,  Vc,  was  determined  experimentally: 

Vc  =  0.40  m/sec 

This  was  done  by  gradually  filling  the  tank  to  the  test  depth  of  0.60  m,  which 
developed  the  equilibrium  bed.  The  current  velocity  was  slowly  increased  until 
weak  movement  was  observed.  A  few  sand  particles  just  began  to  move  in  isolated 
areas  of  the  bed.  The  current  velocity  was  measured  with  the  current  meter  5.0 
cm  above  the  bed.  This  procedure  was  repeated  several  times. 


52 

Comparison  between  the  experimentally  determined  V'c  and  the  values  in  Fig. 
3  reveal  a  close  correspondence.  The  experimental  value  of  V'^  =  0.40  m/sec  is 
slightly  higher  then  that  shown  in  Fig.  3.  This  may  be  due  to: 

(1)  the  precision  of  the  current  meter,  and 

(2)  measuring  the  velocity  at  5.0  cm  above  the  bed. 

Measurements  less  than  5.0  cm  from  the  bed  result  in  the  current  meter 
disturbing  the  bed.  The  V'c  profile  is  shown  along  with  the  velocity  profiles  of 
each  test  condition  in  Fig.  18.  V'c  is  determined  under  current  only  conditions.  A 
hypothetical  tangent  line  from  5.0  cm  above  the  bed  on  the  velocity  profile  curve 
on  Fig.  18  to  the  plane  of  the  bed  (x-axis)  gives  theoretical  values  for  velocity. 
These  values  are  in  the  boundary  layer.  Incipient  motion  occurs  at  velocities  in 
the  boundary  layer.  For  purposes  of  this  study,  the  values  used  are  measurable 
quantities  taken  at  5.0  cm  above  the  bed. 

At  U  =  Vc  one  may  observe  a  slight  movement  in  the  bed.  Close  observation 
revealed  a  few  grains  of  sand  rolling  along  the  bed  in  the  direction  of  the  current. 
At  U  <  V'c  no  sand  moves.  When  the  test  wave  was  run  with  V^  the  measured 
velocity  at  5.0  cm  was  the  same  as  the  measured  velocity  for  current  only.  But 
with  the  wave,  sediment  is  brought  into  suspension,  and  the  sand  bed  mobilizes. 
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At  the  onset  of  each  test,  the  sand  bed  was  leveled  and  the  tank  slowly 
filled  with  water.  Fresh  sand  was  added  to  the  bed  between  tests,  and  manually 
compacted  to  a  thickness  of  0.2  m.  This  was  an  attempt  to  control  the  bed  density 
without  taking  compaction  measurements,  ensuring  the  same  bed  density  for  all 
tests  and  prevented  bed  liquefaction. 

Equihbrium  bed  developed  as  the  tank  filled  with  water.  The  equilibrium 
bed  was  characterized  with  small  sand  ripples,  approximately  2-5  mm  high,  2 
mm  wide  and  about  4-8  cm  apart.  The  bed  remained  in  equilibrium  until  the 
pipes  were  placed  on  the  bed.  Once  the  test  flow  depth  of  60  cm  was  reached,  the 
test  current  was  established  and  measured. 

Pipes  were  placed  in  the  tank  once  the  depth  and  current  were  established. 
For  tests  with  current  and  wave,  the  wave  generator  was  preset  and  started 
immediately  after  the  pipes  were  in  place.  Pipe  placement  methodology  was  the 
same  for  each  test.  Each  pipe  was  placed  in  the  tank  sequentially;  the  first  pipe 
is  the  leading  pipe.  Scour  began  immediately. 

Velocity  profiles  were  measured  and  recorded  for  each  test  at  about  one  meter 
upstream  of  the  leading  pipe.  The  profiles  shown  on  Fig.  18  represent  averages  of 
all  flow  conditions  for  each  appHcable  case.  All  the  profiles  shown  in  Fig.  18  are 
measured  values.  The  wave  gauge  was  located  downstream  of  the  pipe  test  area. 
Fig.  16  shows  the  typical  test  section  layout. 
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Each  test  was  run  between  60  and  90  minutes.  Initially,  scour  occurred  at 
a  high  rate.  Scour  profiles  were  traced  on  sheets  of  acetate  taped  to  the  side 
of  the  tank.  These  scour  profiles  were  taken  periodically,  with  color  coded  lines 
representing  different  times.  Generally,  the  maximum  scour  depth  was  reached 
within  the  first  twenty  minutes.  The  time  profiles  differ  from  the  final  measured 
centerhne  profiles  because  of  side  wall  interference.  These  profiles  indicate  the 
time  when  the  scour  hole  was  at  maximum  depth.  Once  the  maximum  depth 
was  reached,  the  test  was  allowed  to  run  for  twenty  minutes,  ensuring  equilibrium 
scour. 

A  test  run  was  stopped  once  it  became  apparent  that  the  scour  hole  reached 
its  maximum  equilibrium  depth.  The  water  slowly  drained  out  of  the  tank  once 
the  wave  generator  was  shut  down  and  the  flow  control  valve  closed.  The  pipes 
were  carefully  removed  from  the  tank,  without  disturbing  the  scoured  bed.  The 
centerhne  profile  was  measured  with  the  point  gauge  profiler  at  intervals  of  five 
milhmeters  or  one  centimeter.  The  measured  profile  was  plotted  for  each  test. 
The  measured  bed  centerline  profile  and  maximum  scour  depth  for  each  test  is 
included  in  Appendix  II. 
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CHAPTER    VII 
ANALYSIS  OF  RESULTS 


Presentation  of  Data 

The  Maximum  equilibrium  scour  depths  and  relative  scour  depths  which 
were  measured  during  the  tests  in  this  study  are  presented  in  Table  II  and  III, 
respectively.  The  data  show  greater  scour  depth  for  current  when  compared  to 
similar  cases  of  current  and  wave  action,  (oscillatory  motion).  Table  III  shows  the 
relative  scour  depth  Sd/D  decreasing  with  increasing  pipe  diameter.  The  pipe 
spacing  is  the  clear  width  between  outer  pipe  diameters. 

With  two  pipes  in  contact  with  the  bed,  there  were  three  cases  of  no  scour. 
Scour  was  not  evident  when  3.0  cm  pipes  were  exposed  to  current  and  wave  action. 
The  small  pipe  diameter  and  oscillatory  motion  allowed  for  the  pipe  to  capture 
sediment.  The  turbulence  resulting  from  vortex  shedding  was  influenced  when 
an  additional  pipe  was  placed  in  the  wake  of  the  leading  pipe.  Where  no  scour 
occurred,  little  or  no  vortex  shedding  occurred  behind  the  leading  pipe.  There 
was  insufficient  turbulence  in  the  space  between  the  pipes  to  keep  sediment  in 
suspension  for  transport.  Thus,  sediment  was  deposited  between  the  cylinders. 
The  oscillatory  motion  contributes  to  deposition  by  moving  sediment  back  between 
the  pipes.  The  sediment  entrained  between  pipes  settled  out  of  suspension.  In 
currents,  sediment  was  transported  in  the  direction  of  the  current. 
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T^ble  n.  Maximum  Scour  Depth. 


MAXIMUM   SCOUR   DEPTH   -Sd-  (cm)        | 

TEST 
CONFIGURATION 

U  a  0.37  m/t 

U  s  0.26  m/s 

PIPE  DIAMETER  -  D  -  (cm)    | 

3 

6 

9 

3 

6 

9 

—  o 

2.15 

3.75 

4.60 

1.55 

2.20 

2.25 

—  o 

1.45 

3.05 

3.90 

1.20 

2.00 

2.00 

—  OfO 

1.75 

3.25 

4.95 

1.05 

2.15 

2.25 

—  OIO 

1.40 

2.80 

3.70 

no 
scour 

2.05 

2.15 

—  o-o 

2.40 

3.80 

5.15 

1.40 

2,20 

2.25 

—  o-o 

2.15 

3.35 

4.45 

no 
scour 

2.30 

2.30 

—  OfOfO 

1.85 

4.35 

4.65 

1.40 

2.35 

3.25 

—  OfOfO 

1.40 

3.80 

3.60 

0.85 

1.80 

2.25 

^  o»  o  •  o 

2.65 

4.65 

5.10 

1.45 

2.00 

3.75 

^=  0»  O"  0 

1.95 

4.10 

3.95 

1.10 

1.90 

2.50 

—  010 

o 

1.50 

3.05 

4.40 

NOTES  : 

—  0 

currant  only 

o 

0.75 

2.95 

3.95 

p 

2.40 

3.90 

4.35 

=^0 

ctimnt  &  wav« 

—  o»o 
o 

no 
teour 
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Table  m.  Relative  Scour  Depth. 
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In  most  cases  of  multiple  pipes,  the  vortex  shedding  from  the  leading  pipe  is 
influenced  by  the  presence  of  additional  pipes.  This  is  seen  in  the  scour  profiles 
presented  in  Appendix  II.  In  m.ost  cases,  sediment  deposition  occurs  in  the  space 
between  multiple  pipe  configurations.  The  deposition  of  sediment  between  pipes 
is  the  result  of  reduced  turbulence  in  the  wake  of  the  leading  pipe.  Interruption 
of  the  vortex  shedding  by  additional  pipes  results  in  less  turbulence,  or  a  smooth 
wake.  Thus,  sediment  deposition  occurs  in  the  space  between  pipes. 

Subtle  variations  in  scour  profiles  and  relative  scour  depths  are  apparent 
when  comparing  current  only  to  current  and  wave  action  profiles  of  similar  cases. 
Compared  with  current  only,  current  and  wave  action  results  in  a  more  active  bed 
with  less  relative  scour.  This  paradox  is  attributed  to  the  oscillatory  wave  motion, 
disturbing  capability  of  wave  and  wave  energy  transmitted  to  the  bed. 

The  disturbing  capability  of  the  wave  is  described  by  the  wave  steepness, 
[H/gT^).  In  oscillatory  motion,  sediment  is  lifted  into  suspension  and  moved  in 
oscillatory  motion.  The  pipe  spacing  reduces  the  vortex  shedding,  thus  producing 
less  turbulence.  The  wave  pushes  a  particle  back  into  the  scour  hole  where  there 
is  less  fluid  motion.  Sediment  in  suspension  begins  to  settle,  sediment  deposition 
results,  and  the  scour  profile  develops.  The  developing  scour  profile  influences  the 
flow  streamlines. 
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Flow  streamlines  around  multiple  cylinders  during  scour  is  complex.  Stream- 
lines are  in  a  dynamic  state,  influenced  by  the  changing  scour  profile  and  cylinder 
shape.  The  developing  scour  profile  alters  the  flow  streamlines.  The  streamlines 
adjust  to  the  scour  profile.  EquiHbrium  occurs  when  fluid  velocity  becomes  less 
than  the  critical  velocity  of  sediment  movement.  In  equilibrium,  streamlines  are 
uniform  and  the  scour  hole  profile  is  static. 

Data  Verification 

Data  verification  is  achieved  by  comparing  the  single  pipe  case  test  results 
of  this  study  with  results  presented  by  Nalluri  and  Ibrahim  (1986),  as  shown  in 
Figs.  4-7  and  Figs.  19-22.  In  Figs.  19-22  the  single  pipe  restdts  of  this  study 
are  plotted  along  with  the  data  of  Nalluri  and  Ibrahim.  These  plots  indicate  good 
correlation  between  the  single  pipe  cases  of  this  study  and  the  study  of  Nalluri 
and  Ibrahim.  ' 


Figs.  19  and  20  show  the  relative  scour  depth  plotted  against  the  pipe  Froude 
number  and  sediment  Froude  Number,  respectively.  Flow  velocity  is  part  of  the 
Froude  Number,  which  represents  an  inertia-gravity  force  system.  The  relative 
scour  depth  Sd/ D  increases  with  increasing  pipe  and  sediment  Froude  Number. 
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Fig.  21  shows  the  influence  of  the  wave  Ursell  Number,  Ur  =  HL^/y^,  which 
represents  the  relative  magnitude  between  H/L  and  y/L.  Ur  <  I  corresponds 
to  a  region  of  intermediate  to  deep  water.  The  data  from  this  study  represent 
intermediate  to  deep  water  conditions.  The  data  of  Nalluri  and  Ibrahim  (1986) 
are  for  shallow  water  conditions. 

Fig.  22  shows  wave  steepness,  H/gT^,  which  represents  the  disturbing 
capability  of  the  wave.  The  waves  used  in  this  study  are  longer  than  the  waves 
used  by  Nalluri  and  Ibrahim  (1986),  and  therefore  plot  at  lower  Sd/D  values. 

Pipe  Spacing  Relationship 

The  relationship  between  pipe  spacing  and  scour  depth  achieved  with  the 
data  in  this  study  are  dimensionless  plots  of  the  product  of  pipe  Reynolds  Number 
(Nrp)  and  sediment  Froude  Number  {Ns)  against  relative  scour  depth  (Sd/D). 
Scour  is  a  function  of  pipe  geometry,  sediment  size  and  flow  conditions.  The 
Reynolds  Number  and  Froude  Number  represent  apphcable  scaling  laws. 

The  dimensionless  product  of  pipe  Reynolds  Number  and  sediment  Froude 
Number  account  for  the  most  significant  parameters  which  influence  scour.  Plots 
developed  between  Sd/D  and  {Nrp){Ns)  have  a  constant  sediment  Froude 
Number  (Ns)  since  the  sediment  is  the  same  in  all  cases.  The  pipe  Reynolds 
Number  is  variable  and  contributes  to  the  decreasing  Sd/D. 
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Fig.  23-26  are  plots  of  tests  with  U  =  0.37  m/sec.  Fig.  27-30  are  plots 
of  tests  with  U  =  0.26  m/sec.  These  are  plots  of  {Nrp){Ns)  against  (Sd/D)  . 
The  solid  line  in  each  figure  is  the  single  pipe  case,  representing  a  relative  datum 
line.  The  lines  represent  a  linear  least  square  fit  for  a  set  of  three  data  points. 
Statistically  speaking,  more  data  points  are  required  for  definitive  conclusions, 
therefore,  the  lines  represent  only  general  trends. 

Fig.  23  represents  a  trend  where  the  1/2  diameter  pipe  spacing  results  in  less 
scour  than  both  the  single  pipe  case  and  full  diameter  case.  Similar  results  are 
obtained  in  Fig.  24,  current  and  wave.  The  1/2  diameter  pipe  spacing  results  in 
less  scour  than  the  single  pipe  case  and  full  diameter  case.  Comparing  Fig.  23 
with  Fig.  24  shows  a  general  trend  for  greater  relative  scour  in  the  current  only 
case. 

Fig.  25  represents  a  general  trend  where  the  1/2  diameter  pipe  spacing  resvilts 
in  less  relative  scour  than  the  full  diameter  pipe  spacing.  When  compared  to  the 
single  pipe  case,  there  is  no  apparent  trend.  The  full  diameter  spacing  results 
in  greater  scour  than  both  single  pipe  and  1/2  diameter  spacing  cases.  In  Fig. 
26,  the  full  diameter  pipe  spacing  is  about  the  same  as  the  1/2  diameter  spacing. 
Both  the  1/2  and  full  diameter  spacing  result  in  greater  relative  scour  than  the 
single  pipe  case.  Comparing  Fig.  25  with  Fig.  26  shows  a  general  trend  for  greater 
relative  scour  in  the  current  only  case. 
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In  Figs.  23-26,  the  general  trend  is  that  1/2  pipe  diameter  pipe  spacing  results 
in  less  relative  scour  than  the  full  diameter  pipe  spacing.  The  1/2  diameter  pipe 
spacing  for  2-pipe  systems  results  in  less  relative  scour  than  both  the  single  pipe 
case  and  the  full  diameter  case.  The  1/2  diameter  spacing  of  the  3-pipe  system 
appears  to  result  in  greater  relative  scour  than  the  single  pipe  case.  The  full 
diameter  spacing  for  a  3-pipe  system  results  in  greater  relative  scour  than  the 
single  case.  This  may  be  due  to  the  wide  profile  characteristic  of  a  3-pipe  system. 

Figs.  27-30  are  plots  with  U  =  0.26  m/sec.  In  these  cases  the  bed  was 
less  dynamic  than  in  the  set  of  tests  with  U  =  0.37  m/sec.  The  general  trend 
is  similar.  The  1/2  diameter  spacing  results  in  less  relative  scour  than  the  full 
diameter  spacing  for  2-pipe  systems.  In  Fig.  28  the  trend  lines  connect  two  data 
points.  This  is  the  case  where  there  was  no  scour  of  the  2-pipe  system,  3  cm 
diameter  pipes.  The  no  scour  points  are  plotted  at  {Sd/D  —  0.0). 

Figs.  23-30  show  a  trend  for  decreasing  Sd/D  with  increasing  {Nrp){Ns). 
Since  A'5  is  constant,  this  trend  is  due  to  Nrp,  which  represents  increased  inertia. 
Conceptually,  inertia  is  the  tendency  of  a  body  to  resist  acceleration.  Reynolds 
Number  is  an  inertia-viscous  force  ratio.  Increasing  the  pipe  diameters  increases 
the  pipe  Reynolds  Number  (Nrp)  and  results  in  greater  system  inertia. 

In  the  model  test,  pipe  diameters  varied  while  sediment  size  remained 
constant.  The  larger. pipe  diameters  yield  a  higher  pipe  Reynolds  Number.  Table 
II  shows  that  Sd/D  decreases  as  D  increases  for  each  set  of  test  configurations. 
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The  trend  for  decreasing  Sd/ D  as  {Nrp){Ns)  increases  is  attributed  to  increasing 
pipe  diameter,  which  is  a  function  of  Nrp. 


The  scour  profiles  presented  in  Appendix  II  show  a  variation  in  scour  profiles 
obtained  between  the  two  test  currents.  Although  the  3-pipe  system  results  in 
greater  relative  scour  depths,  scour  depth  beneath  additional  pipes  appears  to  be 
less.  This  is  especially  true  for  the  U  =  0.26  m/ sec  cases.  The  general  trend  for 
the  relative  scour  depth  to  be  greater  in  the  3-pipe  system  is  because  the  scour 
hole  profile  is  wider.  Inspection  of  the  3-pipe  configurations  with  U  =  0.26  m/sec 
(Appendix  II)  reveals  many  cases  where  the  third  pipe  behind  the  leading  pipe 
has  minimal  if  any  scour.  When  examining  the  profiles  in  Appendix  II,  note  the 
distortion  of  the  axis  scales. 

The  equilibrium  bed  consisted  of  ripple  formations,  discussed  in  Chapter 
VI.  The  ripples  migrate  in  the  direction  of  the  current.  The  migration  is  not 
visible,  but  apparent  when  comparing  time  profiles.  Tests  with  current  and  wave 
action  resulted  in  a  more  active  bed.  In  these  cases,  sediment  is  transported  from 
upstream  of  the  leading  pipe.  As  the  sediment  is  transported  through  the  pipe  test 
section,  there  is  some  sediment  deposition.  In  some  cases  sediment  was  deposited 
in  the  vicinity  in  front  of  the  leading  pipe.  This  combined  with  the  migration  of 
bed  ripples  accounts  for  sediment  deposition  in  front  of  the  leading  pipe  in  several 
cases.  The  buried  pipes  had  no  apparent  effect  on  stabilizing  the  bed. 
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Scour  Rate 

Initially,  scour  occurs  at  a  high  rate.  Fig.  31  shows  three  time  stages  of 
scour  for  test  number  64.  Test  64  was  conducted  on  a  single  9  cm  pipe  exposed  to 
current  and  wave  action  with  U  =  0.267n/sec.  The  maximum  equiHbrium  scour 
depth  is  established  early,  with  equilibrium  profile  at  a  later  time.  In  most  tests 
the  maximum  equilibrium  scour  depth  was  reached  within  twenty  minutes,  and 
the  equilibrium  scour  profile  within  ninety  minutes. 

Fig.  32  shows  three  time  stages  for  test  number  70.  Test  70  is  a  three 
pipe  system  with  1/2  diameter  spacing  between  pipes.  Test  70  is  exposed  to  the 
same  flow  conditions  as  test  64.  The  maximum  depth  below  the  leading  pipe  is 
established  early.  The  third  pipe  experiences  sediment  deposition.  This  is  due  to 
the  reduced  turbulence  in  the  wake  of  the  second  pipe.  The  rate  of  scour  may  be 
a  function  of  the  effective  sediment  size,  as  well  as  flow  conditions.  Fig.  31  and 
Fig.  32  are  typical  examples  of  the  scour  rates  observed  in  this  study.  Also  note 
that  the  direction  of  flow  in  Fig.  31  and  32  is  from  right  to  left. 
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(a)  Scour  profile  after  10  minutes. 


(b)  Scour  profile  after  40  minutes. 


(c)  Scour  profile  after  90  minutes. 


Fig.  31.  Test  number  64:  scour  profile  developixient. 

{U  =  0.26  m/sec,  L  =  1.0  m,  H  =  6.0  cm,  T  =  1.2  sec) 
(Pipe  diameter  =  9  cm) 
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(a)  Scour  profile  after  10  minutes. 


(b)  Scour  profile  after  40  minutes. 


(c)  Scour  profile  after  90  minutes. 

Fig.  32.  Test  number  70:  scour  profile  development. 

{U  =  0.26  m/sec,  L  =  l.O  m,  H  =  6.0  cm,  T  =  1.2  sec) 
(Pipe  diameter  =  9  cm;  spacing  =  4.5  cm) 
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CHAPTER    VIII 

SUMMARY  OF  CONCLUSIONS 
AND  RECOMMENDATIONS 

One  of  the  major  problems  of  placing  pipelines  in  the  ocean  is  pipeline 
scour,  which  is  a  threat  to  pipeline  stability.  As  more  pipelines  are  placed  in 
the  ocean,  many  adjacent  to  one  another,  scour  will  continue  to  be  a  problem. 
Research  objectives  of  this  study  were  to  determine  the  relationship  between  pipe 
geometry  and  scour  depth,  and  examine  scour  hole  dimensions  under  multiple 
pipes.  Pipeline  scour  was  reviewed,  followed  by  a  theoretical  and  descriptive 
presentation  of  pipeline  scour. 

Seventy-two  model  tests  were  conducted  on  multiple  pipelines  in  contact 
with  a  sand  bed.  Model  tests  were  conducted  in  a  two-dimensional  wave-current 
flume  facility.  Multiple  pipe  systems  were  exposed  to  unidirectional  flow  and 
unidirectional  flow  with  oscillatory  motion.  The  model  testing  investigated  pipe 
spacing  and  its  relationship  to  scour.  A  major  influence  on  scour  was  believed  to 
be  vortex  shedding  in  the  wake  of  the  pipe.  Vortex  shedding  was  disrupted  by 
placing  additional  pipes  behind  the  leading  pipe.  A  scour  profile  was  developed 
for  each  test. 
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Model  testing  resulting  in  the  following  conclusions: 

(1)  For  two  pipe  systems,  1/2  pipe  diameter  pipe  spacing  between  pipes 
results  in  less  scour  depth  than  full  diameter  spacing.  This  is  due  to 
interruption  of  the  vortex  shedding  occurring  behind  the  leading  pipe. 
At  a  full  diameter  spacing,  the  vortex  shedding  has  enough  space  to 
develop.  The  1/2  diameter  spacing  blocks  the  vortex  action,  preventing 
its  full  development. 

(2)  Unidirectional  flow  with  oscillatory  motion  results  in  less  scour  than 
unidirectioneil  flow.  Oscillatory  motion  contributes  to  deposition  by 
moving  sediment  back  into  the  scour  hole,  between  pipes.  Since  there  is 
less  turbulence  between  pipes,  the  sediment  settles  out  of  suspension. 

(3)  The  scour  rate  is  high  initially,  then  decreases  until  equilibrium  is 
reached.  Maximum  scour  depth  was  usually  reached  within  twenty 
minutes  during  the  model  tests,  and  equilibrium  profile  was  attained 
in  about  ninety  minutes. 

(4)  A  relationship  exists  between  relative  scour  depth,  Sd/D,  and  pipe 
Reynolds  Number,  Nrp.  Relative  scour  depth  Sd/D  decreases  with 
increasing  pipe  Reynolds  Number  Nrp. 
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(5)      Wave  characteristics  affected  scour  by  disturbing  and  mobilizing  the  bed. 
The  effect  was  a  mass  transfer  in  the  direction  of  the  current. 


Following  are  recommendations  for  further  research: 

(1)  A  larger  data  base  contributes  in  developing  definitive  conclusions. 
Additional  research  should  concentrate  on  one  set  of  pipe  configurations. 
Collect  additional  data  for  one  set  of  configurations  at  various  pipe 
spacings  for  analysis. 

(2)  Conduct  a  series  of  tests  on  different  sediment  sizes  and  pipe  configura- 
tions. This  may  result  in  an  improved  relationship  between  pipe  spacing 
and  sediment  size,  and  establish  correlation  between  scour  rate  and  sed- 
iment size.  This  will  also  develop  the  relationship  between  the  product 
{Nrp){Ns).  By  varying  pipe  Reynolds  Number  and  sediment  Froude 
Number,  a  relationship  between  relative  scour  depth,  pipe  sizes  and  sed- 
iment sizes  will  be  determined. 

(3)  Investigate  various  waves  and  the  effect  of  the  Keulegan-Carpenter 
number  on  multiple  pipe  systems.  Varying  wave  characteristics  will 
develop  stronger  relationships  on  the  affect  of  waves  on  scour. 
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APPENDIX  I 
SCOUR.FOR  PROGRAM  CODE 


SCOUR. FOR 
This  program  calls  subprogram  WAVE123  which  computes 
the  fluid  velocity  directly  in  front  of  the  pipe. 
It  then  calculates  the  maximum  scour  depth  and  scour 
hole  profile  under  the  pipe.  An  output  file,  PROFILE.DAT, 
is  created.   PROFILE.DAT  is  plotted  using  the  VAX/VMS 
PICSURE  program. 


*  Variables 

*U23   wave  velocity  in  ft/s 

*naterial_constant   dimensionless 

*U_wave     m/s 

*U_current  m/s 

♦U__flow    m/s 

*D_pipe     mm 

♦Gravity   u/s**2 

*K_s  dimensionless 

♦Sand  factor   dimensionless 

*H_inIt    mm 

♦Epsilon    0 

*Number_of_points   Number  of  points  to  be  outputed 

*Number_of_points 

♦Sdmax  mm  Maximum  scour  depth 

*X   mm   Position  along  the  sea  bed  perpendicular  to  the  pipe 

*X_inc  mm  Sampling  increment 

*W    mm   Wave  number 

*H   mm   Distance  between  the  bottom  of  the  pipe  and  the  sea  bed 

*PI        Circumference  constant 

REAL  Material_constant 

&,UZ3 

& ,U_wave 

ft ,U_current 

ft,U_flow 

ft ,D_pipe 

ft, Gravity 

&,K_s 

ft, Sand  factor 

fc,H_inIt 

ft ,Epsilon 

& , Numbe  r_of _points 

ft  rSdmax 

&,X 

ft,x  inc 

&,w~ 

&,H 

ftrPI 
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SCOUR.FOR  PROGRAM  CODE  -  CONTINUED 

PARAMETER    ( PI-3 .14159265 ,Gcavity-9 .81 ) 

CALL   WAVE123(U23) 
WRITE(6,*)'WAVE123    COMPLETED' 

K_s-0.972 

U   wave-UZ3*.3048 


C  Open  file  MAIN. DAT  as  input  to  read  in  the  following  parameters 
OPEN(UNIT-10,riLE-'MAIN.DAT' , STATUS-' UNKNOWN' ) 
READ(10,*)U  current 
READ(10,*)D~pipe 
READ( 10, *)Sand_f actor 

Epsilon"H_init  /  D_pipe 

U^flow  •  U^wave  +  U_current 

Sdmax-K_s*(U_flow**2/Gravity)**0.20*D_pipe**0.800 

W-Sand_factor*Djpipe 

Number_of_points-  30 

OPEN(UNIT-20,riLE-'PROriLE.DAT' ,STATUS-'NEW' ) 

VfRITE(  20,*) 'POSITION      DEPTH' 

X--D_pipe*4 

X_inc-2*ABS(X)/Number_of_points 

X-X-X_inc 

DO   I-l,INT(Number__of_points) 
X-X+X_inc 
IF    (    ABS(X) .LT.(W/2)     )    THEN 

H-Epsilon  "fD  pipe/2  +  Sdmax/2*( l+cos( 2*PI*X/W) ) 

ELSE 

H-Epsilon   •t-D_pipe/2 
ENDIF 

H-H-D   pipe/2 
WRITET6,*)X,H 
WRITE{20,*)X,H 
ENDDO 

CALL    EXIT 
END 
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APPENDIX  II 
TEST  PROFILES 
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APPENDIX  III 
NOTATION 

The  following  symbols  are  used  in  this  paper: 

A  =    projected  area  of  particle  on  plane  perpendicular  to  flow; 

c  =    frictional  drag  coefficient  of  sand  bed; 

Cd  ~    drag  coefficient; 

Ci  =    form  coefficient  defining  effective  surface  area  of  particle; 

D  —    pipe  diameter: 

6.3  —    effective  sediment  diameter; 

(fao  —    mean  particle  diameter; 

Fd  —    drag  force  on  particle; 

g  —    gravitational  acceleration; 

H  =    wave  height; 

h(x)  —    vertical  distance  from  pipe  center  to  bottom  of  scour  hole; 

KC  —    Keulegan-Carpenter  Number; 

A'a  =    sediment  factor; 

L  —    wave  length; 

A'/p  =    pip^  Froude  Number; 

}sTp  —    pipe  Reynolds  number; 

A  .s  —    sediment  Froude  Number; 

i?e  =    Reynolds  Number; 

R.  =    boundary  Rej'nolds  Number; 

Sd  —    maximum  equilibrium  scour  depth; 

5/  —    sand  factor; 

T  =    wave  period; 
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NOTATION  -  CONTINUED 

U  =^  mean  local  fluid  velocity  near  pipe; 

Uc  =  current  velocity; 

I'd  —  local  velocity  at  a  distance  d  above  the  bed; 

l-r  —  wave  Ursell  Number; 

I'w  —  wave  velocity; 

t^ioo  =  fluid  velocity  at  100  cms  above  the  bed; 

U  —  mean  fluid  velocity; 

l\  =  turbulent  shear  velocity; 

\'c  =  critical  velocity  of  sediment  motion; 

w  =  width  of  trench; 

y  =  water  depth; 

Z  =  depth  of  interest; 

A  =  relative  density  of  sediment  in  water; 

e  =  distance  between  pipe  and  bed; 

/i  =  dynamic  viscosity; 

u  =  kinematic  viscosity; 

p  =  fluid  density; 

Pa  =  sediment  density; 

9  —  Shields  parameter; 

T  =  shear  stress; 

T(.  =  critical  shear  stress;  and 

To  =  bed  shear  stress. 
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